Large Force Fluctuations in a Flowing Granular Medium 
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We report the characteristics of the temporal fluctuations in the local force delivered to the wall 
of a 2D hopper by a granular medium flowing through it. The forces are predominantly impulsive 
at all flow rates for which the flow does not permanently jam. The average impulse delivered to 
the wall is much larger than the momentum acquired by a single particle under gravity between 
collisions, reflecting the fact that momentum is transferred to the walls from the bulk of the flow 
by collisions. At values larger than the average impulse, the probability distribution of impulses is 
broad and decays exponentially on the scale of the average impulse, just as it does in static granular 
media. At small impulse values, the probability distribution evolves continuously with flow velocity 
but does not show a clear signature of the transition from purely collisional flow to intermittently 
jamming flows. However, the time interval between collisions tends to a power law distribution, 
P(t) ~ t -3//2 , thus showing a clear dynamical signature of the approach to jamming. 
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Sand, rather than a liquid, is chosen for the contents 
of an hourglass because the rate of efflux from a column 
of sand flowing under gravity does not depend on the 
height to which it is filled. Likewise, in a static column 
of sand, the pressure at a given depth is independent of 
the height of the column above it jy. In both cases, this 
is because the weight of the sand in the interior of the 
column is borne by the lateral containing walls. In the 
case of a static granular column, experiments |2|, ||, Q, 
theory |J ^| , and numerical simulations [§], Q, || have all 
shown that stresses in the bulk are transmitted to the 
walls in a very spatially inhomogeneous fashion. This 
is demonstrated by two independent observations: first, 
the distribution of forces, P{f), at the boundaries of the 
medium is exponential, rather than gaussian j^, [l(J Sec- 
ondly, grains that are highly stressed are organized into 
filamentary, linear structures called 'force chains' that 
carry a large fraction of the stress [||, ||, || . 

In this article we report experiments directed toward 
producing a complementary understanding of a flowing 
granular medium. The force chains in a static medium 
are unstable to perturbations perpendicular to their 
length; when jostled by other grains moving in a flow, 
do they merely become short-lived, or do they melt away 
entirely? Concomitant to this presumed annealing away 
of the stress inhomogeneities, does the broad, exponen- 
tial force distribution in the static case become gaussian 
when grains move and rearrange? The answers to these 
questions have broad consequences for the prospects of 
applying continuum theories to flowing granular matter. 
If indeed the bounding walls of a medium can commu- 
nicate with the interior by transient force chains, then 



'Electronic address: 
t Electronic address: 



aeaswar@smith.edi] 



menon@physics.umass.edr 



any approach to a continuum limit must take into ac- 
count these long length scales. Previous analyses of the 
force distributions in a static granular medium ||, show 
that an exponential at large forces may be obtained from 
scalar models in which the weight of a bead is borne by 
the beads beneath it in some random proportion. While 
force chains do not naturally emerge from these models 
it is tempting to speculate that force chains and broad 
force distributions are related manifestations of the large 
stress inhomogeneities in a granular medium. This point 
of view is also suggested by recent simulations that show 
a simultaneous narrowing of force distributions and a 
blurring of the force chains |sj , as static granular me- 
dia are subjected to compression. On the other hand, a 
recent proposal [ [Ti| for a unified description of jamming 
in thermal as well as non-thermal systems identifies sig- 
natures of the approach to a jammed state in the force 
distribution, P(f). In contrast to earlier models, they 
suggest that loss of mobility is due to the formation of 
force chains, whose presence is most clearly signalled by 
the scarcity of unstressed regions (i.e. a dip or plateau 
in P(f) at / < /, the average force) rather than by the 
exponential tails at large force values. The formation of 
a plateau in P(f) has been also been reported in a recent 
simulation [T^| of grains flowing down an inclined chute. 

In this article we report measurements of the temporal 
fluctuations of the forces delivered to the side walls of 
a hopper (shown in Fig. la) which contains a 2D flow of 
steel balls with diameter d=3. 125mm. The sides of the 
hopper are at a fixed angle of 10° to the vertical, so that 
there are no locations within the cell where the flow stag- 
nates. A transducer (PCB Piezotronics model 209C01) 
sits flush against one of the sidewalls of the hopper and 
measures forces normal to that wall. The diameter of 
the head of the transducer ks d; thus there is only one 
ball against the transducer at any time. For most of the 
measurements we report, the transducer is located at a 



FIG. 2: (a) Impulse histograms on a log-linear scale for 
various V) (all in cm/s): 9.4(X), 11.7(+), 14.7(«), 17.7(4), 
23.0(T), 28.4(A), 33.1(c), 37.6(v), 44.5(A), 50.1(0) and 
60.0(D). The flow velocities correspond to opening sizes rang- 
ing from a — 3d to a — 16d. The data at 9.4 cm/s and 23 
cm/s correspond to the transducer located higher up in the 
flow, (b) Impulse histograms of Fig. 2a scaled to the average 
impulse I for each flow velocity. 
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FIG. 3: (a)Probability distributions, P(t), of the time inter- 
vals r between collisions, on a log-log scale for a flow velocity 
of Vf = 14.7cm/ s corresponding to a = 3.3d. The solid line 
corresponds to a power law P(r) ~ t~' a/2 . The average time 
interval between impulses f , is marked in the figure, (b) P(r) 
on a log- log scale for different flow velocities, Vf (all in cm/s): 
9.4(X), 11.7(+), 14.70), 17.7(*), 23.0(T), 28.4(a), 33.1(c), 
37.6(v), 44.5(A), 50.1(0) and 60.0(D). The corresponding 
opening sizes range from a = 3d to a = 16d. The data at 9.4 
cm/s and 23.0 cm/s are taken with the transducer located 
higher up in the flow. The curves are displaced vertically for 
clarity. The solid line is a power law: P(t) ~ T ~ 3 / 2 
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~ 20msecs, which lies beyond the largest time interval 
for which we are able to obtain good statistics. How- 
ever, if the power-law of r~ 3 / 2 is indeed the asymptotic 
shape of the distribution, then the mean time interval 
t = f P{r)d,T tends to diverge just as in a glass tran- 
sition (even though the average time computed from a 
finite, albeit large, data set shows a relatively innocuous 
dependence on Vf as in Fig. lc). At even slower flow ve- 
locities (the two curves above the solid line in Fig. 3b), 
there are increasingly frequent instances of sustained con- 
tact: at Vf = 11.7cm/ s and 14.7cm/s the percentage of 
the time spent in sustained contact is 9% and 4%, respec- 
tively, compared to 0.3% at Vf — 44.5cm/s. This leads 
to ambiguities in defining the time interval between col- 
lisions possibly resulting in the dashed curves being even 
broader than the power-law of t~ 3 / 2 . A simulation of 
gravity-driven channel flows fll^| also showed a power- 
law distribution of collision times with an exponent of 
close to 3. In that simulation, however, they were able to 
compute the distribution of collision times for particles in 
the bulk, averaged over time, whereas in the experiment 
we only have access to the collisions at a fixed location on 
the wall and cannot keep track of any correlations that 
are convected down the flow with individual balls. 

In summary, our data show that the collisional regime 
extends very close to the threshold of jamming. The 



broad, exponential distribution of forces in the jammed, 
static case extends right through the jamming transition 
and up to the largest flow velocities that we are able 
to explore. The small force end of the force distribu- 
tion does evolve as the flow is slowed down, but the 
indication of an imminent jammed state appears via a 
dynamical quantity, namely, the distribution of collision 
times, P(t). The existence of transient force chains in a 
flowing medium remains unresolved by the measurements 
reported here, though it is clear that if force chains are 
intimately connected with exponential force fluctuations, 
then they must persist in even the most rapid flows we 
have investigated. We are pursuing this issue further by 
direct measurement of spatial correlations between dif- 
ferent points in the flow. 
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